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Abstract 

We carry out a new study of quark mass matrices M u (up-type) and M d (down-type) which are 
Hermitian and have four zero entries, and find a new part of the parameter space which was missed 
in the previous works. We identify two more specific four-zero patterns of M u and M d with fewer 
free parameters, and present two toy flavor-symmetry models which can help realize such special 
and interesting quark flavor structures. We also show that the texture zeros of M u and M d are 
essentially stable against the evolution of energy scales in an analytical way by using the one-loop 
renormalization-group equations. 
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I. INTRODUCTION 


The discovery of the Higgs boson [1] signifies the “completion” of the Standard Model 
(SM) which is not only phenomenologically successful but also theoretically self-consistent. 
In particular, it verifies the Brout-Englert-Higgs mechanism and Yukawa interactions which 
are responsible for the generation of lepton and quark masses. However, the SM is not really 
“complete” in the sense that it cannot explain the origin of neutrino masses, the structures 
of lepton and quark flavors, the asymmetry of matter and antimatter in the Universe, the 
nature of dark matter, etc. Hence one has to go beyond the SM and explore possible new 
physics behind it in order to solve the aforementioned puzzles. 

Here let us focus on the flavor puzzles in the SM. The flavor issues mainly refer to the 
generation of fermion masses, the dynamics of flavor mixing and the origin of CP violation. 
Even within the SM in which all the neutrinos are assumed to be massless, there are thirteen 
free flavor parameters which have to be experimentally determined. On the other hand, one 
is also puzzled by the observed spectra of lepton and quark masses and the observed patterns 
of flavor mixing, which must imply a kind of underlying flavor structure j2J. 

In this paper we restrict ourselves to the flavor issues in the quark sector where there are 
ten free parameters: six quark masses, three flavor mixing angles and one CP-violating phase. 
Thanks to the coexistence of Yukawa interactions and charged-current gauge interactions, 
the flavor and mass bases of three quark families do not coincide with each other, leading 
to the phenomenon of flavor mixing and CP violation. The latter is described by a 3 x 3 
unitary matrix V, the so-called Cabibbo-Kobayashi-Maskawa (CKM) matrix [3], 



( y ud 

V us 

0 

V = 

V cd 

V ca 

V* 


K V td 

V ts 

V tb j 


which can be parameterized in terms of three mixing angles (9 12 , 9 13 , 9 23 ) an d one CP- 
violating phase (5) via the definitions of V us = sin 9 12 cos 0 13 , V ub = sin6 , 13 e -1 ' 5 and V cb = 
cos 9 13 sin # 23 and the unitarity of V itself. As V originates from a mismatch between the 
diagonalizations of the up-type quark mass matrix M u and the down-type one M d , which are 
equivalent to transforming their flavor bases into their mass bases, an attempt to calculate 
the flavor mixing parameters should start from the mass matrices in the flavor basis. In 
view of the experimental results m u <C m c <C rn t , m d <C m s <C m b and 9 13 <C 9 23 <C 9 12 = 


2 




9 C ~ 13°, where 9 C denotes the Cabibbo angle, we believe that the strong hierarchy of three 
flavor mixing angles must be attributed to the strong hierarchy of quark masses. 

Therefore, one is tempted to relate the smallness of three flavor mixing angles with the 
smallness of four independent mass ratios m u /m c , m c /m t , m d /m s and m s /m b . A famous 
relation of this kind is the Gatto-Sartori-Tonin (GST) relation sin 0 12 ~ \f m dl m s PJ- The 


Fritzsch ansatz of quark 

mass matrices [5], 
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( 2 ) 


can easily lead us to the above GST relation. Note that M F and M d possess the parallel 
structures with the same zero entries. Furthermore, they have been taken to be Hermitian 
without loss of generality, since a rotation of the right-handed quark fields does not affect 
any physical results in the SM or its extensions which have no flavor-changing right-handed 
currents. The Fritzsch ansatz totally involves eight independent parameters, and thus it can 
predict two relations among six quark masses and four flavor mixing parameters. However, 
it has been shown that this simple ansatz is in conflict with current experimental data |6). 

One may modify the Fritzsch ansatz by reducing the number of its texture zeros. Given 
a Hermitian or symmetric mass matrix, a pair of its off-diagonal texture zeros are always 
counted as one zero. Hence the Fritzsch ansatz has six nontrivial texture zeros. It has been 
shown that adding nonzero (1,1) or (1,3) entries to M F and M d does not help much [7j, but 
the following Fritzsch-like ansatz is phenomenologically viable mm 
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( 3 ) 


We see that Hermitian M u and M d have the up-down parallelism and four texture zeros. So 
far a lot of interest has been paid to the phenomenological consequences of Eq. (3) [8HTT]. In 
particular, the parameter space of this ansatz was numerically explored in Ref. |12j . where 
a mild hierarchy |5|/A rv_/ B/\B\ 0.24 was found to be favored for both up and down 

sectors. Here we have omitted the subscript “u” and “d” for the relevant parameters, and 
we shall do so again when discussing something common to M u and M d throughout this 
paper. Although there are four complex parameters in Eq. (3), only two linear combinations 
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of the four phases are physical and can simply be denoted as (j) x = arg(C' u ) — arg(C d ) and 
0 2 = arg(R u ) — arg(R d ). It has been found that 0 2 is very close to zero or 2ir [12] . while 
sin0 1 is close to ±1 and its sign can be fixed by r/ u r/ d sin 0, > 0, where the dimensionless 
coefficients r/ u and r] d will be defined in section II. 

In this paper we aim to carry out a new study of the four-zero texture of quark mass 
matrices and improve the previous works in the following aspects: 

• We reexplore the parameter space of M u and M d by taking into account the updated 
values of quark masses and the latest results of the CKM flavor mixing parameters. 
The new analysis leads us to a new part of the parameter space, which is interesting 
but was missed in Ref. H2 and other references. 

• We identify two more specific four-zero patterns of M n and M d with fewer free param¬ 
eters. Namely, there is a kind of parameter correlation in such an ansatz, making the 
exercise of model building much easier. We present two toy flavor-symmetry models 
to realize such special and interesting quark flavor structures. 

• The running behaviors of M u and M d from a superhigh scale down to the electroweak 
scale are studied in an analytical way by using the one-loop renormalization-group 
equations (RGEs), in order to examine whether those texture zeros are stable against 
the evolution of energy scales. We fold that they are essentially stable in the SM. 

The remaining parts of this paper are organized as follows. In section II we first explore the 
complete parameter space of M u and M d and then discuss the relevant phenomenological 
consequences. Particular attention will be paid to some properties of the four-zero texture 
that the previous works did not put emphasis on. Section III is devoted to discussions about 
the special patterns of four-zero quark mass matrices in which some particular relations 
among the finite matrix elements are possible. Two toy flavor-symmetry models, which can 
help realize such interesting patterns, will be presented for the sake of illustration. In section 
IV we derive the one-loop RGE corrections to M u and M d which evolve from a superhigh 
energy scale down to the electroweak scale. Our analytical results show that those texture 
zeros are essentially stable against the evolution of energy scales. As a byproduct, the 
possibility of applying the four-zero texture of quark mass matrices to resolving the strong 
CP problem is also discussed in a brief way. Finally, we summarize our main results and 
make some concluding remarks on the quark flavor issues in section V. 
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II. THE PARAMETER SPACE: RESULTS AND EXPLANATIONS 


Before performing an updated and complete numerical analysis of the parameter space 
of Hermitian M u and M d with four texture zeros, let us briefly reformulate the relations 
between the parameters of M u d and the observable quantities H2. First of all, M ud can be 
transformed into a real symmetric matrix M u d through a phase redefinition: 
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P^MP = 


( 0 

\c\ 

V 0 
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0 ^ 
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( 4 ) 


where the subscript “u” or “d” has been omitted, and P = Diag{l,e with 

< t> B = arg (B) and <f> c = arg(C'). Of course, one may diagonalize M as follows: 


O t MO 
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( 5 ) 


Without loss of generality, we require A and A 3 to be positive. Then \B\, B and \C\ can be 
expressed in terms of A and the three quark mass eigenvalues \ (for 2 = 1,2, 3, corresponding 
to m u ,7n c ,7n t in the up sector or m d , m s ,7n b in the down sector): 


B — X 1 -j- A2 T A 3 — A , 
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In this case the orthogonal matrix O reads 
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( 6 ) 


(7) 


where 7] = ±1, and the emergence of this coefficient can be understood as follows. Since 
A and A 3 have been taken to be positive, A 3 and A 2 must have the opposite signs so as to 
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assure a negative value of the determinant of M , 


Det(M) = -A\C\ 2 = A 1 A 2 A 3 . 


( 8 ) 


When identifying A 123 with the physical quark masses, we use 7] = +1 and —1 to label 
the cases (A 1; A 2 ,A 3 ) = (—m u ,m c ,m t ) and (A 1 ,A 2 ,A 3 ) = (m u , —m c ,m t ) in the np sector, 
respectively. The same labeling is valid for the down sector. 

In terms of quark mass eigenstates, the weak charged-current interactions are written as 


- V, - Im c i) L 7 " V 


AA 

s 

V 6 / 


w* + h.c. , 


(9) 


where the CKM matrix V appears in the form V = O^P*P d O d . The nine elements of V 
can be explicitly expressed as 

Vu, = Ol.Ot + 0 2 " Ote^ + , ( 10 ) 

where (p 1 = (p c — 0 C and (j) 2 = <fr B — (ft Bd , and the subscripts i and a run over (u, c,t ) and 
(d, s, b), respectively. Now it is clear that V depends on four free parameters A u , A d , </> 1 and 
0 2 , after the quark masses are input. With the help of the above analytical results, we are 
able to constrain the parameter space of M u and M d by taking account of the latest values 
of the CKM matrix elements m 

0.22536 ±0.00061 
0.97343 ±0.00015 
0.0405 ±0.0012 

together with the updated values of quark masses at the scale of M z 

m u = 1.38±g;t? MeV , m c = 6380^ MeV , 

m d = 2.82 ± 0.48 MeV , m s = 570 32 MeV , 

In our numerical analysis, we prefer to use |I4J> | V uh \ , \V cb \ and the CP-violating observable 
sin 2/3 as the inputs because their values have been determined to a very good degree of 
accuracy. Here (3 stands for one of the inner angles of the CKM unitarity triangle described 
by the orthogonality relation V* b V ud ± V* b V cd ± V t * b V td = 0 in the complex plane. The three 
inner angles of this triangle are defined as 


|V| = 


/ 


V 


0.97427 ± 0.00014 
0.22522 ± 0.00061 
0.00886 ± 0.00033 


0.00355 ±0.00015 
0.0414 ±0.0012 
0.99914 ±0.00005 


\ 
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m t = 172.1 ± 1.2 GeV 

m b = 286 O 0 go° MeV • 


( 11 ) 


( 12 ) 


a = arg 
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FIG. 1: The allowed regions of A ud , |F? ud |, B ad and |C ud | as constrained by current 
experimental data in the (r/ u , r/ d ) = (+1,+1) case. 


and their experimental values are |13J 

a = (85.4+^)° , sin 2(3 = 0.682 ± 0.019 , 7 = (68.0lf;jj)° . (14) 

Obviously, the uncertainty associated with sin 2 f3 is much smaller than those associated with 
a and 7. The unitarity of V requires a + f3 + 7 = n. 

FIG. 1(a) shows the allowed region of A u and A d , which are rescaled as r u = A n /m t and 
r d = A d /m b , in the (r/ u , r ] d ) = (+ 1 , + 1 ) case. Since the results of r u and r d in the other three 
cases are not quite different from that illustrated in FIG. 1(a), here we just concentrate 
on the (r] u ,r] d ) = (+1,+1) case for the sake of simplicity. Now that r u ~ r d is a quite 
good approximation as shown in FIG. 1(a), we simply use r to denote both r u and r d when 
their difference needs not to be mentioned. We find that the region of r can be roughly 
divided into two parts: (1) r is close to 1 and mainly lies in the range of 0.8 to 0.9; (2) r 
is around 0.5 and mainly ranges from 0.4 to 0.6. These two parts will be referred to as the 
r ~ 1 and r ~ 0.5 regions, respectively, in the following discussions. The reasonableness of 
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FIG. 2: The allowed regions of <p l and (p 2 as constrained by current experimental data in 

the (v u ,Vd) = (i 1 ; ±;L ) cases. 


this treatment will become clear shortly, since the phase parameters (f > 1 and 0 2 behave very 
differently in these two regions. 

The allowed regions of (j) 1 and </> 2 are shown in FIG. [2j where the possibilities of (r/ u , r/ d ) = 
(+1,+1), (+1, —1), (—1,+1) and (-1,-1) have all been considered. Taking the (+1,+1) 
case for example, we find that </> 2 is very close to 2ti and thus its allowed range can also be 
denoted as </> 2 < 0. In comparison, the allowed range of (p 1 is much wider but it can also 
be divided into two parts in a reasonable approximation: 0 1 ~ 0.57T and (p 1 ~ 1.57T. There 
is actually a correlation between r and (f ) l : in the r ~ 1 region (j) 1 ~ 0 .57T holds, and in the 
r ~ 0.5 region (p 1 ~ 1.57T holds. After examining all the four 0 u ,?7 d ) = (±1, ±1) cases, we 
obtain the more general correlation between r and (j) 1 2 as follows: 

r] u r] d sin0 1 > 0 for r ~ 1 ; Vu Vd sin^i <0 for r ~ 0.5 ; 

Vd sin 0 2 < 0 ; Vu Vd cos 0i < 0 • 

Note that only the constraint ?/ d sin0 2 < 0 is numerically exact, and the other three con¬ 
straints serve for good approximations in which most scattered points are satisfied. Such 









TABLE I: The correlation between r and (j) 1 2 in the four (r/ u , r] d ) = (±1, ±1) cases. 



r ~ 1 

r ~ 0.5 

(hu^hd) = (+!.+!) 

(j) 1 ~ 0.5-7T , cos cj) 1 < 0 , cj) 2 < 0 . 

(j) 1 ~ 1.5-7T , COS (j) l < 0 , 4> 2^0. 

(hu^hd) = (+!>-!) 

(j) 1 ~ 1.57T , COS (j) 1 > 0 , (j) 2 > 0 . 

4> 1 ~ 0.57T , COS (j) l > 0 , <^2^0. 

(hu^hd) = (-!>+!) 

(j) 1 ~ 1.57T , COS (j) 1 > 0 , 4> 2 < 0 . 

4> 1 ~ 0.5-7T , COS (j) l > 0 , 4 >2 < 0 . 

(hu^hd) = (-!>-!) 

0! ~ 0.57T , COS ©!<(), 02 > (1 • 

~ 1-57T , COS0! < 0 , 4 >2 > 0 . 


correlative constraints can also be given in a more explicit way, as listed in TABLE I. Finally 
let us point out that the r ~ 1 region and its corresponding parameter correlation found 
here are consistent with the results presented in Ref. [T2], but the r ~ 0.5 region and its 
parameter correlation are our new findings which were missed in the previous works (mainly 
because (p l takes totally different values in this region from our expectation based on its 
values in the r ~ 1 region). 

All the correlative constraints listed in TABLE I can find an explanation once the an¬ 
alytical expression of the CKM matrix V is explicitly presented. No matter whether the 
region r ~ 1 or r ~ 0.5 is concerned, one can easily check that A is close to the mass of the 
third-family quark and thus it is much larger than the masses of the first- and second-family 
quarks. As a result, the orthogonal matrices O u and 0 A can approximate to 



Because m u /m c ~ m c /m t ~ sin 4 9 C and m d /m s ~ m s /m b ~ sin 2 0 C hold, the (1,3) entry 
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of O u is negligibly small but that of 0 d is not. Although the factor m s /(r d m b ) is actually 
much smaller than 1, it is kept in the (3,1) and (3,3) entries of O d since it will play a crucial 
role in explaining the correlation rj d sin <p 2 < 0. 

Given the approximate results of O u and O d in Eq. (16), it is straightforward to calculate 
all the CKM matrix elements by using Eq. (10). We are particularly interested in 



m d 

1 

Vn Vd\ - 


V m s 

1 

V r u (1 r d 

1 r ^~ / 

m d m s f 


V m b TTl b V 


(17) 


m. 


wi„ 


cb 


Among them \V cb \ deserves special attention and can be decomposed as follows: 


\VJ = I ^ (V' b ) - i Im (V' b ) | = ^[Re(E c ;)] 2 + [Im(E c ;)] 2 , 

Re (Kb) - \A U (1 - r d ) ~ ^(1 - r u ) r d ^1 - ^ 1^0 C os0 2 ; ^ 

Im (R C ft) - Jo- - O r d (i - ^ ^0 sin 02, 

where V^ b = has been defined. Clearly, neither |Re (V^) nor Im (VjjJ is allowed to 

be larger than the experimental result \V cb \ ~ 0.04. That is why r u is always nearly equal to 
r d and </> 2 is so close to 0 or 2n. For either r u ~ r d ~ 1 or r u ~ r d ~ 0.5, the fact of cj) 2 ~ 0 
(or 27r) allows us to simplify the expression of |V^ S | to 


IK, 




- 2 i) u i)j 



COS + 



(19) 


It is known that the term y/mjm s itself can £t the experimental value of \V US \ to a good 
degree of accuracy (i.e., the GST relation), and hence one has to control the contribution 
from the smaller term \Jvn u lm c by adjusting the CP-violating phase (p v This observation 
immediately leads to cos0 1 ~ 0, or equivalently (f> 1 ~ 0.57T or 1.57T. As first pointed out in 
Ref. [15], the relation in Eq. (19) is essentially compatible with the orthogonality relation 
V* b KmI+K ftRed +^td = 0 after the latter is rescaled by V* b , leading to the striking prediction 
a ~ (p 1 ~ 0.57T for the corresponding CKM unitarity triangle. Needless to say, this prediction 
is consistent with current experimental data shown in Eq. (14). 
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In order to understand the correlation between the signs of sin d, 2 and those of r/ u d , one 
needs to consider the impact of the CP-violating observable sin 2 (3 on the parameter space 
of M u and M d . Eqs. (10) and (16) allow us to obtain 



Then the definition of /3 in Eq. (13) leads us to 


tan/3 = — 


Re(V cd y; b ) Im(M - Im(^) Re(V td V t * b ) 
MV cd V; b ) R e(V td V t l) + lm(V cd V; b ) Im (V td V t l) 


- Vu Vd sm </q 


m c ™d r 'd m b V 


, , Im m s ] Im(K',) 

X 1 - Vu Vd COS( Pl\ — — I T T 19 


m c 


Given the experimental value of sin2/3 in Eq. (14), we arrive at tan/3 = 0.394 ± 0.015. In 
the r ~ 1 region the first term of Eq. (21) is dominant, and thus r/ u r/ d sin 0 1 is required to be 
positive. Note that this term is at most 0.322, if the values of quark masses in Eq. (12) are 
input. Hence the second term of Eq. (21) has to be positive too. In other words, ?/ d sin</> 2 
should be negative because Im(V^ b ) is proportional to sin</> 2 . Furthermore, r] u r] d cos (f> 1 is 
likely to be negative to enhance the contribution of the second term of Eq. (21) to tan/3. 
When the r ~ 0.5 region is concerned, we find that the second term of Eq. (21) becomes 
important, so 77 d sin</> 2 is still required to be negative. Since this term has a chance to 
saturate the experimental value of tan /3, the first term of Eq. (21) is possible to be negative 
in such a case. In fact, ? 7 u ? 7 d sin <p l must be negative in the r ~ 0.5 region if we take into 
account the constraint from \V ub \. With the help of Eqs. (17) and (18), we have 



+>;„ id s * n 0 1 



( 22 ) 


Taking r d = 0.5 for example, we find that the first term of Eq. (22) is about 0.0044, larger 
than the experimental value \V ub \ ~ 0.0036. In this case the second term of Eq. (22) should 
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FIG. 3: The numerical outputs of V uh \ j | V ch versus sin 2/3 and a versus 7 in the 

(h u? 7 d) = (+ 1 ; + 1 ) case. 


be negative so as to partially offset the contribution from the first term. Therefore, we are 
left with r] u r] d sin^ < 0. However, when the first term is large enough, the second term will 
fail to offset its extra contribution, bringing about a lower bound about 0.3 for r d as shown 
in FIG. 1(a). 

To summarize, we have performed a new numerical analysis of the four-zero ansatz of 
quark mass matrices by using the updated values of quark masses and CKM parameters. 
We find a new part of the parameter space of this ansatz — the r ~ 0.5 region together 
with the relevant correlation between <p 12 and r] ud . We have also explained the salient 
features of the whole parameter space of M u and M d in some analytical approximations. As 
a byproduct, FIG. [3] shows the numerical outputs of V uh |/1 V ch versus sin 2/3 and a versus 
7 in the (r] u ,r] d ) = (+ 1 ,+ 1 ) case. One can see that the uncertainties associated with the 
CP-violating quantities a and 7 remain quite significant, and they mainly originate from 
the uncertainties of m s , r u and r d . I 11 the next section, we shall go back to the quark mass 
matrices themselves to look at their structures and to see whether they assume some special 
patterns with fewer free parameters. 


III. SPECIAL FOUR-ZERO PATTERNS AND MODEL BUILDING 


In the Fritzsch ansatz of quark mass matrices, the first term of Eq. (22) should be replaced 


with (■ m s /m b ) \J™ d /m b , whose size is about 0.00078. This result can easily be understood 
from the trace of M^ d (i.e., A = + X 2 + A 3 ), which gives rise to r d = 1 — 0{m s /m b ) in 


the down sector. Since the upper limit that the second term of Eq. (22) can reach is about 
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0.00236, the experimental value of \V ub \ has no way to be saturated by these two terms in the 
Fritzsch ansatz. When the four-zero texture of M u d is concerned, the existence of nonzero 
(2,2) entries modifies the trace of M ud to the form A + B = A x + A 2 + A 3 . The constraint 
on r d is consequently relaxed, and it actually becomes a free parameter. Given a typical 


value r d = 0.9 in the r ~ 1 region, for example, the first term of Eq. (22) contributes a 
value 0.0014 to \V ub \ such that the experimental result of \V ub \ can be well fitted. In this 
case the magnitude of B is about 0.1A 3 . To avoid a relatively large A 2 from such a large B, 
the three parameters A, B and \B\ must satisfy an approximate geometrical relation p3] up 
to a correction of 0(m 2 /m 3 ): 


\B\ B 


_ _ 

l + O ( —M 

A ~ \B\ 



(23) 

This observation is certainly supported by the numerical results presented in FIG. [lj In 
light of the definition A/m 3 = r, B/m 3 ~ 1 — r holds as a good approximation because of 


|Ai + A 2 | < A 3 = m 3 . Hence Eq. (23) implies \B\/m 3 ~ y/r (1 — r) . In short, the (2,3) 
sectors of M n and M d have the same structure which can be parameterized as 


M <131 ~ A 

u ^ i u 


M ?- 31 


A. 


(24) 


where e ~ ^/(l — r) /r , and its value is about 0.3 in the r ~ 1 region of the parameter 


space. Eq. (24) hints at a common origin of the (2,3) sectors of M u and M d , and thus it can 


be taken as a guideline for model building. However, a numerical analysis shows that such 
an up-down parallelism is slightly broken by the (2,2) entries of quark mass matrices. In 
addition, their (1,2) entries do not share this kind of parallelism, as one can see in FIG. [l](d) . 
With the help of Eq. (j6|, we typically take r u ~ r d ~ 0.9 and illustrate the finite matrix 
elements of M u and M d as follows: 


M ~ A 

± x U - ■ L U 


( 


\ 


0 0.0002 0 

0.0002 0.11 0.31 

0 0.31 1 


\ 


M d ~ A d 


( 


\ 


0 0.005 0 

0.005 0.13 0.31 

0 0.31 1 


\ 


(25) 


It is worth reiterating that the mild hierarchy in the (2,3) sectors of quark mass matrices is 
crucial to fit current experimental data. 

In the r ~ 0.5 region of the parameter space of quark mass matrices M u and M d , there 
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is a particularly interesting case, 


A, = B„ 




(26) 


which deserves special attention. We have verified that these exact equalities are really 
allowed in our numerical calculations. The corresponding parameter space is certainly a 
part of the parameter space restricted by r ~ 0.5. In this special case, the (2,3) sectors of 
M n and M d have a neat form: 


~ A 


t 2m e \ 


mn t 


1 _ ^ 




1 i 

2 m„ 




/ 


2m 0 


m. 


\ 


V 


m. 




A typical numerical illustration of the structures of M u and M d turns out to be 


( 


M — A 

U - u 


0 0.0005 0 

0.0005 1 0.993 

0 0.993 1 


\ 


( 


M d ~ A d 


0 0.013 0 

0.013 1 0.96 

0 0.96 1 


\ 


(27) 


(28) 


One can see that the (2,3) sectors of quark mass matrices are suggestive of an underlying 
flavor symmetry which controls the second and third quark families. 

In fact, the 2 3 permutation symmetry of quark mass matrices, which is quite similar 

to the striking /i <(-)■ r permutation symmetry in the lepton sector |16j . has been conjectured 
long before im Under this simple flavor symmetry the mass matrix takes the form 


M = 


( 0 

c 

\ c 


c 

A 

B 


C^ 

B 

A l 


(29) 


But such a scenario has been ruled out by the present experimental data, as pointed out in 
Ref. [18] . This situation can be easily understood by taking a look at the expression of \V ub \ 
in Eq. where the two terms originate from O u and O d in the following way: 


9 


d 

13 



v A 



(30) 


If there were an exact 2 3 permutation symmetry, both 9\ d and 9f 3 would have to be 

vanishing. However, the second term alone is unable to fit the experimental value of V ub |, 
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as already discussed above. Hence we conclude that quark mass matrices might possess a 
partial 2 -H- 3 permutation symmetry such that 

M 22 = M 33 , M 12 7 ^ M 13 . (31) 

Since there is a large hierarchy between (1,2) and (3,3) entries of M (i.e., M 33 M 12 ), 

the 2 O 3 permutation symmetry can be taken as a starting point for model building, 
and it is broken later on by introducing a small (1,2) entry. Furthermore, the equality 
M 22 = M 23 = M 33 should be a good leading-order approximation. 

Having identified two special patterns of four-zero quark mass matrices, we proceed to 
discuss the model building issues in order to derive them. There are several ways to de¬ 
termine or constrain quark flavor structures, among which flavor symmetries should be the 
most popular and powerful one. So far a number of flavor symmetries, such as the Abelian 
U( 1 ) flavor group [12] and the non-Abelian S(3) flavor group [2D], have been tried in this 
respect. Before introducing a flavor symmetry to realize the above special patterns of quark 
mass matrices, let us discuss what the Hcrmiticity of M u d implies for model building. 

Quark mass matrices originate from the Yukawa interactions and are in general non- 
Hermitian and complex. There are two possibilities of making them Hermitian: (a) a proper 
transformation of the right-handed quark fields, or equivalently a proper choice of the flavor 
basis, as one has done in obtaining Eq. (2) or (3) in the SM or its extensions which have no 
flavor-changing right-handed currents; (b) imposing a reasonable assumption, such as the 
parity symmetry to be discussed soon, on the Lagrangian of Yukawa interactions. Note that 
case (a) is no more favored for our present purpose, because an implementation of possible 
flavor symmetries is also basis-dependent, and hence it is hard to coincide with the chosen 
basis of Hermitian quark mass matrices in most cases. So let us focus on case (b) in the 
following model-building exercises. 

Under the parity symmetry, a flavor theory should be invariant when a left-handed 
fermion field is replaced by its right-handed counterpart (i.e., Vb “^ Vy)> or v i ce versa - 
As for the Yukawa interactions of quark fields, the parity transformation is 

yij^i( H )K + yijK( H )^L *—> + ( 32 ) 

where i and j are the quark flavor indices, and ( H) stands for the vacuum expectation value 
(VEV) of the Higgs field. The invariance of Yukawa interactions under parity transformation 
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TABLE II: The fields relevant for the Yukawa couplings and their charges under I/(1) FN . 


Qi/Qi c 

QHQl 

Ql/Q 3 R 


4>2 

$ 3 /$ 4 

Si 

S 2 

-5 

4 

0 

0 

-1 

1 

-4 

-1 


requires the Yukawa coupling matrix elements to satisfy the condition y i:j = y and hence 
the corresponding quark mass matrix must be Hermitian in the flavor space. We are therefore 
motivated to consider Hermitian quark mass matrices in the framework of the Left-Right 
(LR) symmetric model with an explicit parity symmetry [21] - 

The LR model extends the SM gauge groups to SU(2) L x SU(2) R x U( 1) B _ L , where 
SU(2 ) r is the opposite of SU(2 ) L and acts only on the iso-doublets constituted by the right- 
handed fields, and B—L stands for the baryon number minus the lepton number. All the 
fermion fields are grouped into iso-doublets as follows: 


Qi = 


Uj 


d\ 


Qk = 


u 


R 


d 


T l - 


vi 


L 1 k = 


R 


R 


'R 


(33) 


In the present work we concentrate on the quark sector and leave out the lepton fields L l h 
and L R . At the scale A R which is higher than the electroweak scale, SU(2) R x U{1) B _ L is 
broken to I/(1) Y - The residual SU{2) h and U( 1) Y are exactly the SM gauge groups which 
are subsequently broken by a bi-doublet field <f> under SU(2) h x SU(2) R . 


$ = 


VEV 


<$> = 


H 


(34) 


K, 


The six quarks acquire their masses via their Yukawa interactions with <f>: 

y'ij \ i i i + h.c. . 


(u 


l d l) 


hi 


11 




K 


n 


d{ 


(35) 


R 


In the minimal non-supersymmetric LR model n' has a relative phase as compared with n, 
and this may violate the Hermiticity of quark mass matrices. Hence we prefer to (but not 
necessarily) work in the framework of the supersymmetric (SUSY) LR model [22] • Note 
that the yh term in Eq. (35) will be forbidden by the holography requirement of the 
superpotential in this framework. 

Now that the issue of Hermiticity has been settled, let us continue to build quark mass 
models under certain flavor symmetries in a usual way. We begin with a model that can 
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lead to a four-zero texture of M u and M d in the r ~ 1 region. It is easy to derive the special 
pattern of M ud in Eq. (24) with the help of the Froggatt-Nielsen (FN) mechanism [T^j. The 


point is to introduce a global f/(l) FN symmetry to structure the quark mass matrices. All 
the fields relevant for quark masses and their charges under f/(l) FN are listed in TABLE [IT} 
According to the convention in SUSY, Q l R is represented by its corresponding left-handed 
chiral superfield Q^. In the SUSY LR models, at least two bi-doublets are needed to avoid 
the exact parallelism between M u and M d . In our model four bi-doublets are introduced, 
and their VEVs are written as 


<*i> = I ^ 


($ 3 > = 






Ko 



Key 


v 3 J \ K '4 

In addition, two gauge singlets S l and S 2 are introduced to spontaneously break the I/(1) FN 
flavor symmetry. 

For clarity, let us explore the phenomenological consequences of this model step by step. 
The contribution from <h 1 can be expressed as 


(36) 


V33Qr^iQl + 2/23^R^iQl ^ ^23^R ^ 2/22*3r f ^ 


(37) 


where y 22 and y 33 are real, but y 23 is complex. A is the scale where all the fields associated 
with the FN mechanism reside. The non-renormalizable operators arise from integrating out 
the heavy fields which are not explicitly given in TABLE [TT| and thus they are suppressed 
by A. The key point of the FN mechanism is to assume that the ratios of (S' 1 ) and (S 2 ) to 
A are small quantities which can be generally denoted as e, such that each element of quark 
mass matrices is encoded in a power of e. Here we have identified this small quantity with 


the one in Eq. (24), and thus its magnitude is about 0.3. When S 1 and $ 1 acquire their 
VEVs, the (2,3) sectors of M u and M d are of the form 


(38) 



( V 22 c 2 

2/23 \ 


( 2/22 c 2 

—e \ 

M^’ 3 ^ ~ y 33 n i 

y.3.3 

2/33 

> M d’ 3) ~ 2/33«l 

2/33 

2/33 


2/23 £ 

7 2/33 

1 ) 


V 23 e 

7 2/33 

1 J 


which can reproduce the flavor structure in Eq. (24). There is the exact parallelism between 


up and down quark sectors, because they have the same origin (i.e., from <h 1 here). However, 
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this situation also brings about two phenomenological problems. One of them is that the 
(2,2) entries of M u and M d actually do not respect this exact parallelism, as we have seen 


in Eq. (25). The other problem is that the (2,3) entries of M u and M d should have a phase 


difference, so as to assure 0 2 7 ^ 0 or 27r. 

To address these two problems, let us take account of the contribution from <f > 2 as follows: 

3 3 3,3 3 5,51 


VkQr^QI (jj +!&9r*29l(x) +J/22<7r*2<? 


A 4 


(39) 


This treatment modifies Eq. (38) to the form 


M 


(2,3) 


2/33^1 


/ ^ 22 £ 2 _|_ 1122 K 2 ^4 

1/33 %3 K 1 

Vh c 2(23 ^2 £ 3 

V 2/33 2/33 K 1 

( 2/22 2 1 ^22 K 2 4 

C | / 0 


2/23 ^ _|_ ^23 ^2 £ 3 \ 


2/33 


2/33 K 1 
1 


M, 


(2,3) 


2 / 33^1 


/ 

2/23 c _|_ ^23 ^2 £ 3 \ 


(40) 


2/33 


2/33 K i 


2/33 


2/23 c _|_ 2/23 K> 2 c 3 


2/33 K 1 
1 


/ 


V 2/33 2/33 K i 

If the ratios n l /n 2 and k\/k 2 are close but not exactly equal to each other, the difference 
between the (2,2) entries of M u and M d will be of 0(e 4 ) ~ 0 . 01 , in agreement with the 


numerical result given in Eq. (25). The difference between the (2,3) entries of M u and 
M d seems to be of 0(e 3 ) ~ 0.03 and in conflict with Eq. (25). One may essentially get 


around this problem by assuming that the phase difference between y 23 and y 23 is about 
7 t/ 2, such that the absolute values of M u23 and M d23 only have a negligibly small difference 
of 0(e 5 ) ~ 0.003. But the phase difference between M u23 and M d23 is of 0(e 2 ) ~ 0.1, just 
consistent with the value of (j) 2 illustrated in section II. 

Finally, <h 3 and $4 can offer finite masses for the first quark family through the terms 


y^Q^Ql + y^Q^zQl + y'i2Q^AQl + i&Qfcs&l , 


(41) 


from which we obtain 


^u!2 ~ Mi2I — Vl2 K 3 + 2/l2 K 4 j ^dl2 — Ki2l — 2/l2 K 3 + 2/l2 K 4 ■ 


(42) 


The reason that we arrange $3 and <f > 4 to have the same quantum number is rather simple: 
in this case the phases of M ul2 and M dl2 can be different, such that we are left with a 
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nonzero (f> 1 . In a complete flavor-symmetry model the smallness of M ul2 and M dl2 should 
also be explained via the FN mechanism as we have done for the (2,3) sectors of M u and 
M d . Instead of repeating a similar exercise, here we simply assume that k 3 , k 3 , k a and k / 4 
are much smaller than their counterparts k\ , n 2 and k 2 . Of course, the elements M u 
and M 13 are vanishing as limited by the relevant flavor quantum numbers. 

When it comes to the particular case M 22 = M 33 , a non-Abelian flavor symmetry is needed 
to realize this equality. The simplest candidate of this kind is the S(3) group which has three 
irreducible representations 1, V and 2. The tensor products of these representations can be 
decomposed as follows [23]: 


x, \ / y, \ I XiUo + XnU, 

XI ] = ( X lVl + ^ 2)1 + ( X lV2 - X 2Vl)l' + 

x 2 ) _ \V2 


X lVl X 2V2 


(43) 


X 


~ X 2V 


X y v = I ] , X r x Vi’ = Mi ■ 

V 2 \ x ^/ 2 

The quark fields are organized to be the representations of S( 3) in the following way: 


Q L — 1 > Q R — 1 J 


Ql 


2 , 


Q L ) V Qr 


Q 2 k c 

3c 


2 , 


(44) 


while the bi-doublets introduced and their representations under the S(3) group are: 


$1 ~ 1 , $2 - 1' , 


$3 
C I> 1 


— 2 


(45) 


The VEVs of bi-doublets are specified to be 


<*!> = 


Ki 


« i 


<$ 2 > = 


Kj o 


«2 


($ 3 ) = 


(Co 


($ 4 ) = 


Avo 


In this model the equality of M 22 and M 33 results from the operator 


(46) 




Q 3c 


R 


Qi 


yiiQnMQl + QnMQl] , 


(47) 
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and their values are given by 


-^u22 — Mi.33 — i Md22 ~ ^d33 ~ Vl K l 


(48) 


In comparison, the elements M 23 and M 32 are generated by the operators 


«2 + Qr{*2)Q 



(49) 


% [Qfi( ( ^:i)Qi + Q 3 n c (^)Qi] 


which lead us to 


^u23 — V 2 K 2 + V3 K 3 
^d23 = V2 K 2 + V 3 K 3 


M\i32 ~ 


V2 K 2 4" V3^3 1 
^d32 = ~V2 K 2 + V3 K 3 ' 


(50) 


Notice that the Hermiticity of quark mass matrices as required by the LR symmetry leaves 
us y 2 = — y 2 and r/ 3 = y 3 (i.e., y 2 is imaginary while y 3 is real). If only one of the y 2 and 
y 3 terms exists, 02 will be zero, so both of them are necessary. Another noteworthy point 


is that the operators in Eqs. (47) and (49) are completely independent of each other, so it 


is difficult to understand why M 23 is so close to M 22 and M 33 . We conjecture that these 
two operators are possible to come from the same tensor product in a larger group, so that 
M 22 = M 33 = M 23 can be obtained as the leading-order approximation. 

Finally, let us consider the operators 


V3Qr(^2)Ql + V3Qr(^2)Ql ■ (51) 


They lead us to the nonzero (1,2) entries of quark mass matrices: 



-^u !2 ~ M u2 i — 2 / 3 K 2 , -^d !2 — M(12\ ~ V3 K 2 


(52) 


Note that it is (<1> 3 ) = 0 that ensures the vanishing of M 13 and M 31 . There is also a problem 
that 0 4 equals zero, but it can be overcome by introducing the column vector ($ 4 , <h 5 ) T . 


Similar to (4> 2 , < h 3 ) T , $ 4 acquires its VEV but d> 5 does not. In this case Eq. (52) is modified 
to the form 


M aX2 — .W, l2 | — y 3 K 2 + |/ 3 K 4 , -^d !2 — 1 /<J 21 — V 3 K 2 + V 3 K \ • 


(53) 
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We just need kJk 2 ^ k' a /k' 2 to make (p 1 nonzero. The last remarkable issue is that one 
needs to impose the FN quantum numbers on Ql and Q^, in order to explain why the 
magnitude of M 12 is suppressed by a power of e as compared with those of M 22 and M 23 . 
Such a treatment can also help avoid a large M n arising from the operator I/sQr ($ 1 )<5l- If 
we assign an FN quantum number n to both Ql and Qh, for instance, the contribution of 
this operator will be suppressed by e 2n and thus negligibly small. 

In short, we have identified two special four-zero patterns of quark mass matrices and 
discussed two toy models for realizing them. We should point out that the introduction 
of so many bi-doublet Higgs fields may cause the flavor-changing-neutral-current (FCNC) 
problem. But this problem can be avoided by assuming that the LR symmetry breaks at 
a very high scale and there is just one (two) effective Higgs field(s) (as linear combinations 
of the above Higgs fields) at the low scale in which case we go back to the SM (MSSM) 
situation. Otherwise, we can address this issue by introducing some flavon fields located 
at a superhigh energy scale to play the role of bi-doublets as multiple representations of 
the flavor symmetries. In this case, we do not need Higgs fields other than the usual ones 
which have already been required for other purposes rather than the flavor physics. After 
integrating out the flavon fields, there will be no trace of the flavor physics except that the 
Yukawa couplings have been constrained by the flavor symmetries. This way of preventing 
the flavor physics from disturbing the other physics is widely used in flavor-symmetry models 
for the lepton sector |24j . 


IV. ON THE STABILITY OF THE FOUR-ZERO TEXTURE 

As shown in section III, the four-zero texture of quark mass matrices may result from 
an underlying flavor symmetry. But the failure in discovering any new physics of this kind 
indicates that it is likely to reside in a superhigh energy scale, such as the grand unification 
theory (GUT) scale. This means that a flavor-symmetry model should be built somewhere 
far above the electroweak scale and the RGE running effects have to be taken into account 
when its phenomenological consequences are confronted with the experimental data at low 
energies [23]. One may follow two equivalent ways to consider the evolution of energy scales, 
provided there is no new physics between the flavor symmetry scale A FS and the electroweak 
scale M z [26] : (a) the first step is to figure out quark masses and flavor mixing parameters 
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from M u and M d at A FS , and the second step is to run these physical quantities down to 
M z via their RGEs; (b) the first step is to evolve M u and M d from A FS down to M z via 
their RGEs, and the second step is to calculate quark masses and flavor mixing parameters 
from the corresponding quark mass matrices at M z . Here we take advantage of way (b) to 
examine the stability of texture zeros of M u and M d against the evolution of energy scales 
in an analytical way. The RGE effect on the Fritzsch texture of quark mass matrices has 
been studied in a similar way [26, [27]. 

At the one-loop level, the RGEs of the quark Yukawa coupling matrices in the SM can 
be written as 


16tt 2 


d Y q (t) 
d t 


- G q (i)l + T(t) 1 


Y*(t) . 


(54) 


where t = ln(/i/M z ), and the subscript “q” stands for “u” and “d”. The contributions of the 


charged leptons and neutrinos to Eq. (54) have been omitted, because they are negligibly 
small in the SM. Denoting the VEV of the Higgs field as v, we can express the four-zero 
texture of Y u and Y d at A FS as follows: 


( 0 c u 0 ^ 


Ti(^-Fs) — ^Mu(A fs ) — 


c u b u b u 


V 


0 b„ a 


( 0 c d 0 ^ 


> 7d(Aps) — ~dT d (A FS ) — 


»/ 


Cd h b d 


V 


0 b* d a d 


(55) 


Without loss of generality for CP violation, we have chosen b u and c u to be real in Eq. (55). 
The terms G q (t) and T(t) read 


G u — G d + g\ — 8gl + -g% + — g\ , 


T = 3Tr(i;n t + >d>d t 


(56) 


which arise from quantum corrections to the quark and Higgs field strengths, respectively. 
They are flavor-blind, and thus proportional to the identity matrix in the flavor space. Since 
their effects are simply to rescale quark mass matrices as a whole at a lower energy scale, 
they will be dropped for the moment. Namely, we are mainly concerned about the first 


term in Eq. (54): S u — —S A — Y u Yj — Y d Y d , which governs the nonlinear evolution of Y 


Defining H Q = Y q Y T, let us rewrite Eq. (54) by dropping its G q (t) and T{t ) terms: 


32tt 2 d H q 
3 d t 


S q H q + H q S q 


(57) 
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In a good approximation S' can be expressed as 


/ 


q — _ q ~ 

°u — °d — 


0 0 

0 A 2 A 3 


\ 


( 58 ) 


^0 Ag A 1 y 

where A x = A 2 = and A 3 = b u (a u + 6 U ). Then we solve the differential 


equations in Eq. (57) and obtain 


( 


H U (M Z ) ~ 


2 J Pi 

- Cu&uP “ a 7 Ta^ c ’ 


: u a u (a u 6 u - 6 2 


, P -1 

e u fc uP+ a7Ta; c ’ 


■ u b u (a u b u - 


\ 


- A „2 o P P (aX - b 2 


c u + A 2 p 2 — 2 


Ai+A 2 


U^U 


A 3 p 2 


2 a ^ ^ (a u & u -& 2 )' 


A lP 2 - 2 


, (59) 


V ” A, + A 2 v “ “ v 

where the elements denoted as “• • • ” can be directly read off by considering the Hermiticity 
of H u , and p describes the RGE running effects from A FS to M z : 


p = exp 


32vr 2 


y 2 (t')dt' 


(60) 


U FS 


Here t FS = In (A FS /M Z ), and y{t') is the Yukawa coupling eigenvalue of the top quark which 
evolves according to 

= {l y * - G «) yl • (61) 

For illustration, p ~ 0.9 when A FS ~ 10 15 GeV, as shown in FIG. Hi On the other hand, 


/ 


H d (M z ) ~ 


9 " 1 

12 . ,2 


P" 1 -! 


Ai + A 2 


c d b dP 1 - a , A ~ c d ( A i & d - A 3 b d 
12 \ 2 , P ~ P 


P -i — 1 

c d ft dP _1 - A , A c d ( a A - A 3 b d ; 


A ! + A 2 


(A 2 &d ^ A 3 & d ) 


or+ \b d r+b 2 p - 2 + 


o-2 


Ai + A 2 


A 3 ^a d + b^j (b d + b d ) - 2A 1 ( j & d | 2 + 6; 


( a d + & d) P~ 1+ A 1 + A 2 

X ( a d + 2 |k d | 2 + ^d) A 3 


A. 


( a d + l & dl 
( a d + &d) ( b d + b d ) 


a 3 
-2 ^-1 


2 V“2 1 ^ ~ P 


A d + A 2 


2 y^d ^ 
( 62 ) 


l^dl 2 
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FIG. 4: An illustration of the changes of <7 and p with the scale A FS in the SM. 


The RGE-corrected quark mass matrices can then be extracted from Eqs. (59) and (62): 


/ 


M u( M z) ~ 


l P - 1 

b,,p -x-« 

Ai + A 2 


u ( a J>n - b l 


\ 


b - p + a7+xG u (“ u * u - 

(“A-*") 


M d( M z) - &d v 


where 


c d b d P x + 


Cd 


1 


1° 6 >" + yAy( A A- A *«) 


b d p 1 + r , a 1 < A 3«d - 


a dP 


Ai + A 2 

P - 1 -! 
Ai + A 2 


(Aa&d ~ A 2 a d ) 


(j a = exp 


167T 2 


[3 y 2 (t') - G q (t')] d t' 


(63) 


(64) 


is the overall rescaling factor of quark mass matrices brought back from the G q (t) and T(t) 


terms of Eq. (54) that were tentatively dropped in Eq. (57). Apparently, M d (M z ) is not 
Hermitian any more, because the RGE of F d (f) does not respect Hermiticity. To illustrate, 
the numerical changes of er u and cr d with the scale Apg are shown in FIG. [4] in the framework 
of the SM. Of course, the above analytical results can exactly reproduce those obtained in 
Ref. [26] for the Fritzsch ansatz of quark mass matrices when b u and b d are switched off. 
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Note that the geometrical relation in Eq. (23) can be reexpressed as (a q 6 q — 6 q )/(A 1 + 


A 2 ) ~ m 2 /m 3 . Hence in each entry of the (2,3) sector of M U (M Z ) the second term is 
suppressed by a factor proportional to (1 — p) mjm t < 10~ 3 as compared with the Erst 
term. As for M d (M z ). let us take its (2,3) entry as an example to look at the corresponding 
RGE correction. Because of the parallelism between (a u , b u , b u ) and (a d , Re(6 d ), b d ), we find 

Ao a,. 


A 3 a d - AiRe(6 d ) = A 1 Re(b d ) 
- AiRe(fc d ) 


_Ai Re(b d ) 
'A3 Ou _ i 

A K . 


= Re(6 d ) ( ajj u - b 


(65) 


So the real part of the second term of M d23 at M z is suppressed by a factor proportional 
to (1 — p) m c /m t < 10~ 3 as compared with the real part of its first term. In other words, 
the real part of M d23 is approximately equal to p _1 Re(& d ) at M z . According to our phase 


assignment in Eq. (55), 0 2 = arg(6 u ) — arg(6 d ) = — arg(6 d ) holds. Hence the phase of b d is 
equal to — </> 2 and must be close to 0 or —2n. In the r ~ 1 region where A 2 is much smaller 
than A x , it is easy to see that the imaginary part of the (2,3) entry of M d (M z ) is about 
Im(5 d ). That means arg(6 d ) ~ Im(6 d )/Re(6 d ) is rescaled by p due to the RGE effects, or 
equivalently 

02(Afs) — P l< t > 2(Mz) • (66) 

In a word, the four texture zeros of quark mass matrices are essentially stable against the 
evolution of energy scales. To be more specific, M u and M d develop the overall factors a u 
and <7 d during their running from A FS down to M Zl respectively; and their hnite entries 
(a u ,6 u ,fe u ) and (a d , Re(5 d ), b d ) are rescaled by p and p _1 , respectively. 

To illustrate the RGE-induced corrections, let us give a numerical example to compare 


between Eq. (55) at A FS and Eq. (63) at M z . We first figure out the values of quark masses 
and flavor mixing parameters at A FS ~ 10 11 GeV by solving the one-loop RGEs numerically: 


rn u = 0.69 MeV , 
m d = 1.4 MeV , 
K,\ = 0.225 , 


m c = 320 MeV , 
m s = 29.1 MeV , 
IV^I = 0.0458 , 


m t = 95.6 GeV ; 
m b = 1.3 GeV ; 
IVJ = 0.00387 , 


(67) 


and the value of sin 2/3 is almost unchanged from M z to A FS (or vice versa) within the 
accuracy that we need. The choice of this specific scale is for two simple reasons: on the 
one hand, it is expected to be around the canonical seesaw f28j and leptogenesis |29] scales; 
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on the other hand, it is close to the energy scale relevant for the possible vacuum stability 
issue of the SM [3D]. Therefore, 


/ 


^u( A Fs) — 10 


-1 


0 9 x 1(T 4 0 

0.6 1.7 

4.9 


\ 


/ 


^d( A Fs) — 10 


-3 


V. 

0 0.04 e _1 ‘ 67i 0 

1.0 2.4 e 014i 

6.7 


V" 


\ 


( 68 ) 


In comparison, the corresponding quark mass matrices at the electroweak scale are 

/ 0 


M u( M z) - 10 <? u v 


V- 

/ 0 


9 x 10 -4 0 N 

0.6p - 8 • 10 -3 (p - 1) 1.7p + 3 x 10" 3 (p - 1) 

4.9p - 10" 3 (p - 1) J 


M d (M z ) ~ 10 a d v 


0.04 e _1 ' 67i 

„-i 


0 


\ 


l.Op" 1 - (0.2 + O.li) £ 2.4p _1 + 0.3i - 0.04e 

\ ■ ■ ■ (2.4 - 0.3i) p" 1 + (0.05 + 0.04i) e 6.7p“ 4 + (0.01 + O.li) e J 


(69) 


where e = p -1 — 1 is a small value of (9(0.1) or much smaller. This numerical exercise 
confirms our qualitative analysis made above. In particular, the imaginary part of the (2,3) 
entry of M d (M z ) is really independent of p, and its real part is proportional to p _1 . 

Now let us turn to the running behaviors of quark masses and flavor mixing parameters. 
Since c u is negligibly small in magnitude as compared with a u , b u and b u , the invariants of 
the (2,3) submatrix of M U (A FS ) and M U (M Z ) lead us to 

™ c ( a fs) + "h( A Fs) - v (a u + 6 U ) , 

™ c ( a fs) "h( A Fs) - v 2 ( a uK ~ bl) ; 

^ C ( M z) + ™ t ( M z) - a n v («u + K) P , 

m c (M z ) m t {M z ) ~ a 2 u v 2 (a u b u - b 2 ^ p . (70) 

These relations indicate that m c and m t change with the energy scale in the following way: 


"h( A Fs) - o'u P m t( M z) , ™ c ( a fs) - ™ c ( M z) 


(71) 
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When c u is concerned, a similar trick yields m u ( A FS ) ~ a u 1 m u (M z ). It is easy to verify that 
the similar relations hold in the down sector: 


™&( a fs ) ~ a d P m b( M z ) , " 0 ( a fs ) - °d m s( M z ) , 


(72) 


and m d ( A FS ) ~ a d l m d {M z ). These results clearly show that the mass ratios m u /m c and 
m d/ m s are essentially free from the RGE corrections. 

To see how the flavor mixing parameters evolve from A FS down to M z , we take a new 


look at Eq. (17). Above all, the dimensionless parameters r u and r d are independent of the 
energy scale to a good degree of accuracy. The reason is simply that rn t (or m b ) and A u (or 
A d ) have nearly the same running behaviors, as one can see from Eqs. (63), (71) and (72). It 
is also straightforward to conclude that V us is stable against the evolution of energy scales. 
In view of r u ~ r d and 0 2 ~ 0, we arrive at the approximation 


iKbl - yjix - r n) r 6 


1 hd "A 

2 r d m b 


l sm i 


(73) 


Given Eqs. ( 66 ) and (72), the running behavior of \V cb \ turns out to be 


IF*(a F s)I — p |r cb (M z ) . 

With the help of this result and Eq. we immediately obtain 

|V u6 (A fs )| = p-'\V ub (M z )\ . 


(74) 


(75) 


In addition, Eq. (21) tells us that f3 is nearly scale-independent. It is easy to check that a 


and 7 , the other two inner angles of the CKM unitarity triangle, are also free from the RGE 
corrections at the one-loop level [51] . 

The above results can simply be translated into the ones for three flavor mixing angles 
and one CP-violating phase in the standard parametrization of the CKM matrix: 


^12( A Fs) — ^12 {Mz) > $23( A Fs) — P 1 ^23 i^z) ) 

$i3( a fs) — P 1 ^ l i3 {Mz) > ^( a fs) — 8{M Z ) ■ 


(76) 


Of course, a, f3 and 7 are all the functions of 5 in this parametrization. As for the Jarl- 
skog invariant J = cos 0 12 sin 0 12 cos 2 0 13 sin 0 13 cos 0 23 sin 0 23 sin 5 [32], it is easy to arrive at 
J7(A fs ) ~ p~ 2 J(M z ) in the same approximation. Such a rephasing-invariant measure of 
weak CP violation is actually tiny, only about 3 x 10 -5 at M z . 
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Finally, let us briefly comment on a possible implication of the loss of Hermiticity of M d 
running from A PS down to M z . We conjecture that it might have something to do with 
the strong CP problem [33], which is put forward due to the unnatural smallness of the 
parameter 9 = dq CD + 0q FD . Here 0q CD is H ie coefficient of the CP-violating term in the 
QCD Lagrangian [33], 

C„ = ; (77) 

and Oqpd comes from the quark flavor sector, 

#qfd = arg [Det (M u M d )] . (78) 

The experimental upper bound of 9 is at the 10 -11 level [33] . in sharp contrast with a natural 
value of 0(1) from a theoretical point of view. The demand for explaining why 9 is so tiny 
poses the strong CP problem. An attractive solution for this problem is the Peccei-Quinn 
mechanism [36] in which an anomalous U( 1) symmetry is introduced to ensure a complete 
cancellation between ^qcd and #q F d- Another competitive strategy is to remove 9qcd by 
imposing a spontaneously broken P or CP symmetry (e.g., in the LR symmetric model), 
and to keep the second term vanishing in the meantime m Eg. 

Being Hermitian, the four-zero texture of quark mass matrices automatically satisfies the 
requirement arg [Det ( M u M d )\ = 0 at a superhigh energy scale A FS . Nevertheless, the RGE 


effects can render M d (M z ) non-Hermitian as shown in Eq. (63). Since the strong CP term 
begins to take effect at the scale of about 260 MeV where the QCD vacuum transforms, 
nonzero arg [Det (M d )] at or below M z will contribute to 9 in spite of arg [Det (M d )] = 0 at 
A fs . Given the explicit form of M U (M Z ) and M d (M z ) in Eq. (63), one may calculate its 
contribution to 9 as follows: 

•'- 1 - 1 A 3 Im(6 d ) 


0q FD = arg [DetM d (M z )] ~ arctan 


P 


i-i 


(1 ~r) (1 -p) sin 0 2 . (79) 


Ai + A 2 «d P~ 

Although 02 i s very small, it cannot be exactly zero as shown in our numerical analysis. 
Given r = 0.9 and A FS = 1 TeV, for instance, Eq. (79) leads us to a value of 0( 10~ 5 ), much 
larger than the upper bound of 9. One way out of this problem is to fine-tune the value of 
r. But the possibility of r ~ 1 has phenomenologically been ruled out, as discussed at the 
beginning of section III. If the parallelism between the forms of M n and M d is given up, the 
situation will change. For example, in a flavor basis with M u being diagonal, the value of 
60FD was estimated to be of 0( 10" 16 ) in Ref. [39]. In short, it seems difficult to directly 







employ the four-zero texture of quark mass matrices to solve the strong CP problem in the 
scenario of spontaneous CP violation. But a more detailed study of this issue is needed 
before a firm conclusion can be achieved. 

V. SUMMARY 

We have carried out a new study of the four-zero texture of Hermitian quark mass matrices 
and improved the previous works in several aspects. In our numerical analysis what really 
matters is that we have found a new part of the parameter space, corresponding to A ~ \B\ ~ 
B (or r ~ 0.5), and confirmed the known part corresponding to A > \B\ > B (or r ~ 1). 
In particular, the exact equality between A and B is allowed, and this opens an interesting 
window for model building. We want to emphasize that the newly found parameter space 
is phenomenologically different from the already known: since the former allows the (near) 
equality of mass matrix entries — a characteristic of non-Abelian flavor symmetries which 
are very popular in the lepton sector EH, it provides a possibility of unifying the description 
of quarks and leptons with the same flavor symmetries and this will be discussed elsewhere. 

We have identified two special four-zero patterns of quark mass matrices and constructed 
two toy flavor-symmetry models to realize them. One of the patterns possesses a mild 
hierarchy A ~ e\B\ ~ e 2 B with e being about 0.3, and it can be obtained with the help of 
the FN mechanism. The other pattern assumes A = B, which can be realized by means of the 
S(3) flavor symmetry. Both of them show a similarity between the (2,3) sectors of M u and 
M d , indicating that the latter could have the same origin. We have done two model-building 
exercises in the SUSY LR framework with an explicit parity symmetry, which ensures the 
Hermiticity of quark mass matrices at the flavor symmetry scale A FS . 

We have also studied the RGE effects on the four-zero texture of quark mass matrices 
in an analytical way, from A FS down to the electroweak scale M z . Our results show that 
the texture zeros of M u and M d are essentially stable against the evolution of energy scales, 
but their finite entries are rescaled due to the RGE-induced corrections. An interesting 
consequence of the RGE running is the loss of the Hermiticity of M d at M z in the SM. 
As a byproduct, the possibility of applying the four-zero texture of quark mass matrices to 
resolving the strong CP problem has been discussed in a very brief way. 

Although the predictive power of texture zeros has recently been questioned in the lepton 
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sector [3D] , they remain useful in the quark sector to understand the correlation between the 
hierarchy of quark masses and that of flavor mixing angles. We remark that possible flavor 
symmetries are behind possible texture zeros, and they are phenomenologically important 
to probe the underlying flavor structure before a complete flavor theory is developed. 

Note added. While our paper was being finished, we noticed a new preprint [31] in which 
a systematic survey of possible texture zeros of quark mass matrices was done but the 
four-zero texture of Hermitian quark mass matrices with the up-down parallelism was not 
explicitly discussed. 
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